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SUMMARY  (U) 

A  software  package  for  onboard  acoustic  data-processing  and  the  statistical 
analysis  of  array  beam-noise  has  been  implemented  on  an  IBM  3033  computer 
system.  The  package  both  assesses  data  quality  and  measures  the  noise-field's 
statistical  characteristics  and  directionality.  Many  of  the  outputs  are  sufficiently 
general  to  have  application  to  other  typesof  data  or  to  satisfy  other  objectives.  The 
various  outputs  are  described  and  illustrated  by  results  from  measurements  of  the 
ambient  noise  by  a  towed  array.  A  simplified  run  procedureisincluded  to  facilitate 
the  use  of  the  software  package  by  the  analyst. 
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1 .  INTRODUCTION 

Over  a  period  of  several  years,  experience  gained  through  participation  in 
measurement  exercises  and  the  resulting  data  analyses  has  led  to  the 
development  of  a  "standard"  software  package  for  onboard  statistical  analysis 
of  time-series  data.  Although  in  the  past,  the  data  have  been  of  beam-noise 
acquired  by  towed  arrays,  many  of  the  products  are  sufficiently  general  to 
apply  to  data  acquired  by  other  systems,  even  a  single  hydrophone.  The 
purpose  of  this  report  is  to  describe  the  capabilities  and  products  of  this 
software  package.  This  package  is  called  the  Towed  Array  Noise  Analysis 
Package  (TANAP) . 

The  software  package  was  developed  with  two  goals  in  mind.  The  first  was  to 
enable  an  assessment  of  measurement  system  performance  and  data  quality.  This 
requires  routines  that  can  address  the  performance  of  the  sonar  system  as  a 
measurement  tool  and  others  that  can  assess  the  character  of  the  ambient-noise 
field  at  the  time  of  the  measurement.  In  other  words,  how  did  the  condition 
of  the  measurement  tool  affect  the  measurements,  and  was  the  noise-field 
sufficiently  "well  behaved"  to  enable  meaningful  results  to  be  obtained?  The 
second  objective  was  to  provide  products  that  could  be  useful  in  the 
characterization  and  the  analysis  of  the  noise-field.  Some  products  satisfy 
both  purposes  but  are  discussed  only  once.  Some,  such  as  the  mean  level,  are 
commonly  used  and  need  little  or  no  explanation,  so  that  the  discussions  of 
them  are  appropriately  brief.  For  less  well-known  products,  the  discussions 
are  sufficient  to  give  only  an  idea  of  their  form,  utilization,  and  utility; 
more  thorough  discussions  can  be  obtained  from  the  references. 

Section  2  of  this  report  discusses  the  input  data  requirements  of  this 
software  package.  Section  3  discusses  the  statistics  and  products  generated 
that  are  useful  for  assessing  measurement  system  performance  and  data  quality. 
Section  4  discusses  the  statistics  and  products  generated  that  are  useful  to 
analyse  and  describe  the  noise-field  and  the  beam-noise  data. 

Finally,  Appendix  I  contains  a  brief  description  of  the  programs  and  sub¬ 
routines  in  the  beam-noise  processing  software  and  Appendix  II  gives  an 
example  of  a  simplified  run  procedure. 


2.  PRELIMINARY  DATA  PROCESSING  AND  OVERVIEW 

The  experimental  procedure  by  which  the  data  are  collected  and  analysed  is  as 
follows.  Long  duration  (2  to  12  h)  tows  are  made  on  one  heading  to  acquire 
data  to  determine  the  statistics  of  towed-array  beam-noise.  Data  are  also 
collected  while  on  various  array  headings  (called  a  noise  polygon)  to  obtain 
estimates  of  noise-field  spatial  statistics,  ie ,  noise-field  horizontal 
directionality. 

The  data  from  each  hydrophone  are  time  sampled  and  processed  by  a  Fast  Fourier 
Transform  (FFT)  to  get  time-series  of  frequency  spectra  for  each  hydrophone. 
The  spectral  data  for  all  hydrophones  for  each  time  increment  are  then 
processed  by  an  FFT  across  the  array  (with  zero  filling  to  be  a  power  of  2)  to 
get  beam-noise  in  each  spectral  bin.  This  procedure  is  illustrated  in 
figure  1  for  an  array  of  40  hydrophones.  The  resulting  spectrum  analysed 
beam-noise  time-3eries  of  complex  Fourier  beam-noise  coefficients  are  stored 
on  disk  and  are  the  input  data  to  the  beam-noise  processing  and  analysis 
software. 

The  data-processing  system  has  two  purposes.  The  first  purpose  is  to  produce 
information  (called  "quality-assessment  products",  and  described  in  Section  3) 
with  which  to  assess  the  acoustic  performance  of  the  measurement  system  and  to 
determine  whether  the  acoustic  field  is  sufficiently  well-behaved  to  vield 
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Figure  1.  Data  processing  and  reduction 


COMPUTER 


Figure  2.  The  data-acqulsltlon ,  processing,  and  analysis  system 
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meaningful  results,  ie,  to  assess  the  quality  of  the  data.  The  second 
purpose  is  to  produce  and  display  statistical  parameters  (called  "analysis 
products",  and  described  in  Section  4)  that  can  provide  estimates  of  the 
temporal  and  spatial  statistics  of  the  noise-field  (or,  more  accurately,  of 
the  influence  of  the  noise-field  on  the  acoustical  output  of  the  array) .  Some 
of  the  data-processing  products  serve  for  both  assessment  and  analysis. 

Table  1  lists  the  data-processing  products  in  the  order  of  their  appearance. 
The  table  shows  whether  they  are  used  primarily  for  quality  assessment  or  for 
analysis,  or  for  both.  It  also  lists  the  figures  used  as  examples  in 
Chapters  3  and  4. 

The  first  step  in  the  data-processing  (see  figure  2)  is  to  apply  the 
calibration  corrections.  This  produces  ambient-noise  levels  in  dB//uPaJ /Hz . 
An  inverse  FFT  is  then  applied  to  the  calibrated  beam-noise  data  to  provide 
calibrated  hydrophone  data.  When  a  strong  signal  is  present  these  data  can  be 
used  to  calculate  the  phase  angle  across  the  array  relative  to  a  common 
hydrophone  and  normalized  for  the  direction  of  the  source  and  plotted  versus 
hydrophone  number  (see  figure  4). 

Statistical  calculations  for  each  time-series  are  then  undertaken  (see 
figure  5  in  Section  3  below) .  These  calculations  are  performed  on  all  beams 
including  the  virtual  beams  (ie,  those  formed  by  the  FFT  beamformers  that 
correspond  to  phase  shifts  greater  than  that  for  an  endfire  beam)  and  on  two 
hydrophones:  the  one  with  the  median  noise  level  and  another  selected  for 
comparison.  The  results  of  these  statistical  calculations  provide  the  rest  of 
the  data-processing  products  (see  figures  3,  5,  6,  and  7  in  Section  3  and 
figures  8,  9,  10,  11,  12,  and  13  in  Section  4).  The  ways  in  which  these 
products  are  used  for  assessing  the  data  quality  or  for  analysing  the 
noise-field  are  discussed  below  in  Sections  3  and  4  respectivelv. 


3.  ASSESSMENT  OF  TOWED  ARRAY  SYSTEM'S  PERFORMANCE  AND  DATA  QUALITY 

The  measured  data  are  the  basic  foundation  on  which  rests  the  interpretation 
of  results  and  formulation  of  conclusions.  When  the  quality  of  the  data  is 
low  or  uncertain,  the  level  of  confidence  in  the  correct  interpretation  of  the 
results  and  formulation  of  conclusions  must  also  be  low.  Data  quality  can  be 
a  severe  problem  when  the  phenomenon  being  measured  is  as  highly  variable  as 
ambient  noise.  The  severity  is  compounded  when  the  measurement  system  is  a 
towed  line  array,  which  always  has  its  own  interfering  source  (towshipl  and  is 
not  constrained  to  remain  linear,  horizontal,  or  on  a  given  heading. 

The  current  ability  of  the  data-processing  system  to  evaluate  and  filter  out 
self-noise  is  assessed  before  starting  the  measurements  at  a  site.  This 
procedure  is  described  below  in  Section  3.1.  After  the  measurements  have  been 
obtained  and  before  they  are  analysed,  their  quality  is  evaluated.  This 
procedure  is  described  below  in  Section  3.2.  Both  procedures  make  use  of  the 
data-processing  products  listed  in  Table  1. 

3.1  Products 

The  products  and  statistics  which  are  used  to  assess  the  measurement 
system's  performance  and  the  quality  of  the  data  are  listed  below: 
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(a)  To  indicate  contaminating  sources  or  poor  beamforming 

Listings  of  the  average  noise-power  levels,  the  geometric  mean  intensity 
levels,  and  the  differences  between  them.  Single  listings  for  each 
analysis  frequency  in  each  time-period,  as  shown  in  the  AVGPR,  AVO  and 
PRDIF  columns  in  figure  5,  and  plotted  as  in  figures  3(a)  and  3(b)  for 
hydrophones  and  beams  respectively. 

(b)  To  indicate  an  erratic  beamformer  or  hydrophone  output 

Standard  deviations  of  all  beams  and  two  hydrophones.  Single  listing 
for  each  analysis  frequency  in  each  time-period,  as  shown  in  the  STDEV 
column  in  figure  5  and  as  plotted  in  figure  3.  The  beam  numbers  on  the 
outsides  of  the  vertical  dashed  lines  correspond  to  virtual  beams.  The 
utility  of  the  virtual  beams  for  system  performance  and  data  quality 
assessment  is  discussed  later  in  this  section. 

(c)  To  indicate  faulty  hydrophone  channels  or  irratic  beamforming 

The  median  and  average  phase  angle  for  all  hydrophones  normalized  to  the 
steering  direction  of  a  dominant  source  and  the  standard  deviation  of 
phase  angle,  as  shown  in  figure  5  and  plotted  in  figure  4. 

(d)  To  test  for  randomness  in  each  data  set 

Results  of  the  number-of-runs  test  for  all  beams  and  two  hydrophones. 
Single  listing  for  each  analysis  frequency  in  each  time-period,  as  shown 
in  the  ZRUN  column  in  figure  5. 

(e)  To  test  for  transients  and  "outliers"  in  the  data  set 

Results  of  the  mean-square  successive  difference  test  for  all  beams  and 
two  hydrophones.  Single  listing  for  each  analysis  frequency  in  each 

time-period,  as  shown  in  the  ZMSSD  column  of  figure  5. 

(f)  To  detect  linear  trends 

Kendall's  rank  correlation  coefficients  for  all  beams  and 
two  hydrophones.  Single  listing  for  each  analysis  frequency  in  each 

time-period,  as  shown  in  the  ZTAU  column  of  figure  5. 

(g)  To  test  whether  the  towship  noise  is  dominating  the  sidelobes 

Spearman's  rank  correlation  matrix.  Single  matrix  for  each  analysis 
frequency  in  each  time-period,  as  shown  in  figure  6. 

(h)  To  spot  persistent  anomalous  patterns  in  the  noise-level  response 

Polar  plots  of  the  50  percentile  (median)  noise  levels  on  all  beams 
receiving  acoustic  energy.  Done  for  each  analysis  frequency  for  each 
time-period  of  the  directionality  polygons,  as  shown  in  figure  7.  The 
median  level  for  each  beam-noise  time-period  is  also  listed  in  the 
MEDIAN  column  of  figure  5. 


Additional  discussion  of  how  these  statistics  can  be  used  to  assess  the 
performance  of  the  measurement  system  and  data  quality  are  given  in 
references  l,  2  and  3. 
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TABLE  1.  PRODUCTS  OF  THE  DATA- PROCESSING  SYSTEM 


Quality 

assessment 


Description 


Plots  of  fl)  average  noise-power  levels;  (2)  geometric  mean  Incensltv  levels. 3'  median 
level.,  (a)  standard  deviation,  of  the  level.  as  a  function  of  all  bvdrcohones.  Done  ‘or  - 
each  analysis  frequency  in  each  time-period. 

Plots  of  (1)  average  noise-power  levels;  (2)  geometric  mean  intenaltv  levels;  'ji  median 
levels,  U>  standard  deviations  of  the  level,  as  a  function  of  all  beama.  Done  ‘or  each 
analysis  fr.au.ncy  In  ..ch  cln.-p.rlod,  Dlff.rmc,  b.tv..n  (11  and  (2)  ar.  soi-.tln,, 
plotted. 

Hoca  of  (1)  average  and  n.dlan  phas.  and  standard  d.vt.rlon  of  th.  ph.».  a,  ,  'unction  o< 
hydrophone  for  a  high  level  source. 

["•'‘i*”’-  'l  and  90  p.rc.ntlla.  of  nets,  levala  at  all  baama  and  tvo  hydrophones.  ; 
single  listing  for  each  aralyala  frequencv  in  each  time-period. 

Caoparrlc  can  mtanafty  l.v.la  <d.clb.l  av.ra.e.l  on  all  bsui  and  rvo  hydrophones. 

Single  listing  for  each  analysis  frequencv  tn  each  time-period. 

Average  noise-power  levels  on  all  beams  And  two  hydrophones.  Single  listing  for  each 
analysis  fraquency  In  each  time-period. 

Standard  deviations  of  all  beams  and  two  hydrophones.  Single  listing  for  each  aralvals 
fraquency  in  each  time-period. 

Skew  and  kurtoals  for  all  beama  and  tvo  hydrophones.  Single  listing  for  each  ..nalvsis 
'requency  In  each  tine-period.  , 

Results  of  the  number-of-mns  test  for  all  beams  and  tvo  hydrophones.  Single  listing 
for  each  analysis  frequency  in  each  time-period. 

e!!U}f*wMh#  “an~>(»uar*  successive  difference  test  for  all  beams  and  tvo  hydrophones. 

Single  listing  for  each  analyst*  frequency  tn  each  ttma-^aTlod.  | 

Rends  11* a  rank  correlation  coefficients  for  all  heams  and  tvo  hydrophones.  Single  j 

listing  for  each  analysis  frequency  In  each  time-period.  I 

Difference.  b.tw.ao  ch.  .v.r.g.  nola.-pov.r  level  .nd  th.  «, metric  ... n  Int.n.ltv 
l.v.l  on  .11  b..u  .nd  evo  hydrophon...  Slngl.  U.tlng  for  ,.ch  .n.lysl,  fr.qn.ncv 
In  each  tlma-period. 

Spearmen ' .rank  corr.l.tlou  -trig.  Slngl.  matrix  for  ..ch  an.ly.l,  fr.qu.ncv  In  ,«ch 
c  ua*— period . 

Plots  Of  the  cumulative  distribution  functions  of  noise  levels  for  selected  beams  or 
groups  of  baama.  Done  for  each  analysis  frequency  for  each  time-period  and  for  the 
combined  data  sat  from  each  directionality  polygon. 

Polar  plots  of  th*  50  percentila  (median)  noise  levels  on  all  baama  receiving  acoustic 
energy.  Dona  for  each  analyala  frequency  for  eech  time-period  of  the  directions lltv 
polygons. 

Plot,  of  high  r.golutlon  sp.ti.1  ap.ctr.  obt.ln.d  frn»  th.  b«.»-nol».  dst.  bv  th.  V5: 
high  resolution  algorithm. 

Polar  plots  of  the  estimates  of  the  horizontal  directionality  derived  from  the  50  percentile 
(median)  notes  levels  by  deconvolving  tha  beam  response  and  resolving  the  ambiguities.  Dons 
for  each  analysis  frequency  for  the  combined  data  eat  from  each  directionality  polygon. 

Polar  plots  of  tha  estimates  of  horizontal  directionality  derivad  from  the  10,  25,  50 
(median)  percentile  noise  levala  for  beams.  Done  for  each  analysis  freouenev  for  the 
ccnblned  data  ast  from  each  directionality  polygon. 

Oanldlr.cHon.l  noli,  l.v.l.,  Slngl.  vain,  for  ..ch  .n.lr.t.  fr.qu.ncv,  obtained  fron  the 
combined  data  eat  from  each  directionality  polygon. 

Plot,  of  th.  a.lnuch.l  enlaotropy  cumulation  dl.trlhutlon  function  (MCOr) .  Don.  for  each 
an.ly.l.  fr.qu.ncy  for  th.  at.tl.tlc.  tor.  .nd  for  the  conbln.d  data  »et  fron  ..ch 
directionality  program. 

P‘°t"  °(  **‘"  ‘®Ptov«.nt  v.rnia  .rr.y  heading  for  •  .p.clfl.d  detection  sector, 

called  array  heading  roses. 
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TOP  CURVES : 

AVCPR  C 

REDIAR - C 

avg . c 

FOR  ABSOLUTE  LEVELS.  ROD  *XV 

BOTTOM  CURVES' 

stdev - Q 


09  TO  ABOVE 


TOP  CURVES >  _ 

AVGPR  ® 

MEOIAR - ® 

AVC . ® 

FOR  ABSOLUTE  LEVELS.  ROO  'tr09  TO  ABOVE 
BOTTOM  CURVES'  _ 


10  20  30  40  SO  60  14  37  50 


PHONE  NUMBER  BERM  NUMBER 


X  =  4S.EQ  O0S 


FREQ.  =  300-0 

X  :  33. SS  DBS 


(1)  average  noise-power  levels, 

(2)  geometric  mean  intensity  levels, 

(3)  median  levels, 

(4)  standard  deviations  of  the  levels.  (Levels  are  in  dB//uPa1 2 3 4/Hz. )  As 
functions  of  all  hydrophones  or  o'  all  beams  for  each  analysis 
frequency. 

(a)  Phone  level  plots  (b)  Beam  level  plots 

Figure  3.  A  data-processing  output 
3.2  Additional  description  of  statistical  tests 

Several  products  and  statistics  for  the  assessment  of  the  measurement 
system's  performance  and  data  quality  are  listed  in  Table  1,  most  of  which 
require  no  further  explanation.  However,  the  following  four  tests  are  not 
as  universally  well  known  and  will  be  described  in  more  detail  below: 

•  Numbers  of  runs  test 

•  Mean-square  successive  difference  test 
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•  Kendall's  rank  correlation  coefficient 


•  Spearman's  rank  correlation  matrix. 


PHONE  NUMBER 


Figure  4.  A  data-processing  output  :  median,  average,  standard  deviation  and 
percent  of  deviation  in  relative  hydrophone  phase  angle  across  the 
array 


3.2.1  Number-of-runs  test 

The  number-of-runs  in  any  time-period  of  beam-noise  measurements  is  a 
non-parametric  statistic  (ie,  a  statistic  that  is  independent  of  the 
actual  distribution  function  for  the  data  set  under  investigation).  The 
number-of-runs  is  used  as  one  measure  of  randomness  in  a  series  of 
observations.  The  attribute  under  consideration  is  randomness  in  the 
arrangement  of  beam-noise  levels  in  a  time-period  and  not  in  the 
distribution  of  beam-noise  level  magnitudes.  Consequently,  the 
number-of-runs  statistic  indicates  the  temporal  randomness  of  the 
process  that  generated  the  noise  levels  observed  at  the  beamformer 
output  (ie,  the  randomness  of  the  beam-noise  data  in  the  time  domain). 
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KEY  TO  FIGURE  5 


rhe  heading  I'StS  tne  fQl towing 


OATE 
TACTIC 
ST  TIME 
FREQUENCY 
SOUND  SPEEO 
LO*G.  LAT 
sample  size 
LEG 

HEADING 
GAIN 
SI0EL08E 
DEPTHS 
WEIGHT 
ARRAY 
«  AVG 


Date 

Code  t or  the  exercise  site 
atari  time  ot  series 
Analysis  frequency 
In  met  res/ second 

Longitude,  Latitude  of  beginning  of  series 

Number  of  data  in  time-series 

Leg  of  direction*) ity  polygon 

Aosolute  heading  of  array 

Gain  of  signa I -condi tioning  unit 

Assumed  sideioOe  level  of  deaas 

Depths  of  the  array  (■).  Depth  or  bottoa  (a  or  fa? 

Saapl ing  window 

Low- frequency  or  high-frequency  array 
Nuafcer  of  eleaents  averaged  -  T  iae  *  bandwidth 
product  (T*8W) 


Lines  i  to  64  refer  to  beaas.  Lines  6S  and  66  refer  to  hydrophones, 
line  65  being  the  tiae  series  of  aedian  level  data  selected  froa  all 
hydrophones  and  line  66  being  a  single  hydrophone  selected  for 
comparison.  The  coluans  refer  to  the  statistics  calculated  over  the 
saapled  tiae  series,  as  follows. 


The  columns  list,  froa  left  to  right. 


BEAM 

RhOG 

THOG 

THDG 

BW 

Beaa  or  hydrophone  nuaoer 

Heading  (")  of  ambiguous  beaa  pairs  relative  to 
array  orientation 

True  heading  (°)  of  aabiguous  beam  pairs  (right) 

True  heading  (°)  of  aabiguous  beaa  pairs  (left) 
Beaawidth  at  3  dB 

10* 

10  percentile  noise  level  (dB) 

25* 

25  percentile  noise  level  (dB) 

MEOIAN 

50  percentile  (MeJian)  noise  level  (dB) 

75* 

75  percenti le  noise  level  (dB) 

90* 

90  percentile  noise  level  (dB) 

AVG 

Geoaetric  aean  intensity  level  (dB  average) 

AVGPR 

Average  noise-power  level  (dB) 

STOEV 

Standard  deviation  (dB) 

SKEW 

Skew  (dB) 

KURT 

Kurtosis  (dfl) 

dB 

Nuaber  of  beaa  data-saaples  used  to  calculate 
statistics 

ZRUM 

Nuabar  of  runs  test 

ZMSSD 

Mean-square  successive  difference  test 

ZTAU 

Kendall's  rank  correlation  coefficient 

PRO  IF 

Power  difference  (AVGPR  -  AVG)  (dfl) 

The  text  at  foot  reads  as  follows: 

for  beaas  Q  to  □  the  aean  value  of  the  aedian  levels  is  □  dB,  the 
standard  deviation  is  □  dB 

For  beaas  a  to  a  the  aean  STOEV  is  o  dB,  the  standard  deviation  of 
STDEV  is  □  dB. 

For  all  beaas  the  dyneaic  range  of  the  aedian  levels  is  □  d8. 

The  oaniJirectional  noise  level  estiaated  froa  the  aedian  levels  of 
beaas  □  to  □  is  □  dB. 

The  nuaber  of  virtual  beaas  used  is  a  ,  the  virtual  beaa  power  level 
is  a  dB. 

Beaa  a  is  the  aedian  phone. 

Beaa  □  is  phone  □ 

Processed  on  ooo 
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Figure  6.  A  data-processing  output:  Spearman’s  rank  correlation  matrix 
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Figure  7.  A  data-processing  output:  polar  beam  level  plots 
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In  order  to  actually  calculate  the  number-of-runs  statistic (ref . 4) , 
beam-noise  levels  in  the  time-series  must  be  transformed  into 
three  categories,  according  to  whether  they  are: 

■  above  the  median  beam-noise  level  (A's) 

•  below  the  median  beam-noise  level  (B's) 

•  equal  to  the  median  beam-noise  level  (M's). 

Individual  beam-noise  levels  in  the  time-series  are  then  replaced  by  the 
symbols  A,  B  or  M.  Subsequently,  samples  equal  to  the  median  value 
(M's)  are  eliminated  from  the  data  set.  Any  unbroken  sequence  of  like 
kinds  (eg  AA,  BBB,  etc.)  in  the  revised  time-series  is  counted  as  one 
run.  The  number-of-runs  in  the  time-series,  regardless  of  the  length  of 
the  run,  is  then  counted.  The  algorithm  relates  the  number-of-runs 
actually  counted  (RUN)  to  the  distribution  of  runs  expected  for  any 
randomly  distributed  time-series  containing  (N  )  A  elements  and  (N0) 

B  elements.  Assuming  that  all  the  elements  of  the  time-series  have  the 
same  probability  of  occupying  any  position  in  the  temporal  sequence,  a 
probability  distribution  for  each  realizable  permutation  of  the  A's  and 
B's  can  be  determined.  Mathematically,  the  expected  distribution  of  the 
number-of-runs  in  a  large  sample  drawn  from  a  randomly-distributed 
population  can  be  approximated  by  a  normal  distribution  of  the  variable 
R,  with  the  mean  and  standard  deviation  values  given  by: 


R 


2nanb 

NA  +  *B 


+  1 


2NaV2V’b  - 

("a  +  V*  (na 


na- 


+  N_ 


-  1) 


The  algorithm  employed  in  the  software  calculates  the  quantity  ZRUN  (see 
figure  5,  column  17),  which  expresses  the  number-of-runs  actually 
counted  (RUN)  in  terms  of  the  statistics  for  the  corresponding  normal 
distribution: 


ZRUN 


That  is,  ZRUN  is  the  deviation  of  the  actual  number-of-runs  from  the 
expected  mean,  normalized  in  terms  of  the  expected  standard  deviation. 
Beam-noise  time-series  with  too  few  or  too  many  runs  (ie,  those  with 
large  values  of  the  parameter  ZRUN)  are  suspected  to  be  products  of 
non-random  processes.  It  should  be  noted  that  a  completely 
deterministic  signal  -  such  as  an  acoustic  tone  -  could  produce  too  few 
or  too  many  runs,  depending  on  its  frequency,  since  the  sampling 
interval  is  generally  constant. 

The  ZRUN  term  can  be  used  directly  to  eliminate  beam-noise  data  sets 
that  appear  to  have  too  few  or  too  many  runs.  For  example,  one  could 
arbitrarily  state  that  any  time-series  with  ZRUN  >  2  should  be 
eliminated  from  further  consideration.  However,  since  the  expected 
distribution  of  the  number-of-runs  is  normal,  approximately  4.5Z  of  the 
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valid  data  would  be  rejected  by  such  a  criterion.  An  alternate  and 
perhaps  more  appropriate  screening  method  would  eliminate  data  sets 
based  on  the  confidence  level  associated  with  the  hypothesis  that  the 
measured  number-of-runs  (RUN)  is,  in  fact,  a  member  of  the  normal 
distribution  (R,  0*) . 

3.2.2  Mean-Square  Successive  Difference  Test 

The  Mean-Square  Successive  Difference  (MSSD)  Test  attempts  to  verify 
that  each  measurement  of  beam-noise  level  in  any  particular  time-series 
is  independent  of  the  next  (successive)  measurement  in  that  time-series. 
Unlike  the  Number-of-Runs  Test,  the  MSSD  Test  (figure  5,  Column  18)  uses 
the  numerical  value  of  each  measured  noise  level  as  well  as  its 

relationship  to  the  adjacent  (successive)  measurement  in  the 
time-series.  This  test  is  useful  in  identifying  contamination  of 

time-series  data  due  to  transients  of  both  acoustic  and  non-acoustic 

origin. 

The  Mean-Square  Successive  Difference  Test  is  based  on  the  assumption 
that  the  measurements  in  any  time-series  are  drawn  from  a  normally 
distributed  population  with  unknown  mean  and  variance.  Under  this 
assumption,  two  methods  can  be  used  to  estimate  the  population's 

variance,  a2.  One  is  the  unbiased  estimator: 


i-N 


Lj 

i-t 


The  other  is: 


i-N-1 


i-1 


one-half  of  the  mean  of  the  successive  differences  squared.  If  the 
time-series  data  are  dependent  or  are  not  drawn  from  a  normally 
distributed  population,  the  two  calculations  should  produce  different 
estimates  of  the  variance.  For  example,  if  the  beam-noise  time-series 
measurements  are  accumulated  during  periods  containing  steadily 
increasing  or  decreasing  noise  levels  (ie,  an  upward  or  downward 
sloping  ramp),  the  differences  between  consecutive  samples  would  be 
smaller  than  those  obtained  for  a  completely  random  noise-field.  Hence 
the  d2/2  estimate  of  the  variance  would  be  smaller  than  that  calculated 
using  the  unbiased  estimator. 

When  the  beam-noise  measurements  contain  strong,  constant  signals,  the 
two  estimates  may  or  may  not  be  significantly  different.  If  the  strong 
signal  persists  throughout  the  entire  measurement  interval,  the  two 
estimates  would  be  nearly  identical,  since  all  the  levels  are  nearlv 

uniform.  However,  if  the  strong  signal  is  present  only  during  a 
significant  fraction  of  the  measurement  interval,  the  results  produced 
by  the  two  variance  estimators  will  be  substantially  different.  For 
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these  reasons,  the  MSSD  Test  may  be  viewed  as  a  detection  device  for 
certain  kinds  of  transient  signals  in  the  time-series  data. 

The  algorithm  employed  calculates  the  test  statistic  ZMSSD,  where: 


ZMSSD 


d2/2s2  -  1 
[  (N-2)  /  (N2— 1 )  ]** 


If  the  assumptions  concerning  the  beam-noise  measurements  are  valid,  the 
expected  distribution  of  the  test  statistic  ZMSSD  can  be  approximated  by 
a  normal  distribution  with  a  mean  value  of  zero  and  a  unit  variance. 

The  exact  distribution  of  this  test  statistic  has  been  tabulated  for 
uncorrelated  observations  with  sizes  ranging  from  4  to  60  samples 
(ref. 2).  For  sample  sizes  greater  than  six,  the  normal  distribution  is 
a  very  good  approximation  at  the  0.05  and  the  0.95  significance  levels. 
However,  the  normal  distribution  is  not  a  good  approximation  of  the 
exact  distribution  at  the  0.001,  0.01,  0.99,  and  0.999  significance 
levels  until  the  sample  size  is  greater  than  60.  Nonetheless,  the 
normal  distribution  is  a  useful  approximation,  since  it  always  yields 
conservative  estimates.  To  indicate  the  magnitude  of  the  differences, 
critical  values  of  ZMSSD  were  calculated  at  various  confidence  levels 
using  the  normal  approximation.  Table  2  compares  the  results  with  the 
exact  values  (from  reference  5)  for  a  sample  size  of  60: 


TABLE  2.  COMPARISON  OF  ZMSSD  RESULTS  WITH  EXACT  VALUES 


Significance  levelB 


Critical  values 
based  on 

0.001 

0.95 

0.99 

0.999 

Approximate  normal 
distribution 

-3.091 

-2.327 

-1.645 

1.645 

2.327 

3.091 

Exact  distribution 

-3.013 

-2.306 

-1.649 

1.649 

2.306 

3.014 

Since  the  accuracy  of  the  normal  distribution  approximation  improves  as 
the  sample  size  increases,  the  approximation  should  be  more  than 
adequate  for  beam-noise  sample  sizes  of  60  or  more. 

MSSD  Test  results  (ie,  the  quantity  ZMSSD)  are  interpreted  and  used  in 
much  the  same  manner  as  the  Number-of-Runs  Test  results  (ZRUN).  The 
quantity  ZMSSD  can  be  used  directly  to  eliminate  data  sets  from  further 
processing  actions  or  it  may  be  used  to  test  the  hypothesis  that  the 
beam-level  measurements  contained  in  any  time-series  are  distributed  in 
a  random  manner  (ie,  are  not  correlated).  In  the  latter  case, 
tabulated  critical  values  for  the  appropriate  normal  distribution  are 
required  to  confirm  (or  reject)  the  statistical  hypothesis. 

3.2.3  Kendall’s  rank  correlation  coefficient 

The  Kendall  rank  correlation  coefficient  is  a  non-parametric 
(distribution  free)  statistic  that  provides  a  measure  of  the  correlation 
between  two  se  of  ranked  observations (ref . 6) .  Kendalls'  rank 
correlation  coef i.cient  is  used  to  measure  the  correlation  between  each 
time-series  of  beam-noise  measurements  and  a  set  of  raorotonically 
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increasing  integers  (1,  2,  3 . N) .  The  monotonically  increasing  data 

set  used  for  each  correlation  assessment  is  identical  to  the  ranked 
order  of  any  continuously  increasing  function.  Therefore,  increasing 
trends  in  any  time-series  of  beam-noise  measurements  result  in  positive 
values  of  the  correlation  coefficient,  while  decreasing  trends  produce 
negative  values.  The  absolute  magnitude  of  the  correlation  coefficient 
provides  an  indication  of  the  duration  of  the  trending  portion  of  the 
time-series  relative  to  the  length  of  the  sample  under  consideration. 
Although  the  concept  is  similar  in  some  respects  to  the  number-of-runs 
statistic,  the  Kendall  rank  correlation  coefficient  uses  much  more  of 
the  information  present  in  the  time-series  and  is  considerably  more 
powerful  (in  a  mathematical  sense)  than  the  ZRUN  and  ZMSSD  results 
described  previously  when  testing  for  a  linear  trend. 

In  order  to  calculate  Kendall's  rank  correlation  coefficient,  each 
observation  in  the  time-series  must  be  replaced  by  its  ranking  relative 
to  all  other  measurements  in  that  particular  time-series.  As  used  here, 
the  term  'rank'  implies  that  all  measurements  in  the  time-series  can  be 
rearranged  in  a  monotonically  increasing  sequence  from  the  minimum  to 
the  maximum  value  (ie,  ordered  according  to  the  magnitude  of  the 
measurements).  Each  measurement  value  is  then  assigned  a  numerical  rank 
and,  in  turn,  the  numerical  ranks  are  substituted  into  the  original 
time-series  in  place  of  the  measured  values.  Multiple  occurrence  of 
certain  values  in  the  original  time-series  (ie,  tied  data  points) 
complicates  the  process,  since  it  requires  somewhat  more  bookkeeping, 
but  the  end  result  is  essentially  the  same.  Ties  are  treated  in  the 
usual  manner  and  are  assigned  the  same  fractional  rank.  For  example,  if 
there  is  an  n-way  tie  for  the  rank  r^,  the  numerical  rank  assigned  to 
each  of  these  values  is: 

ri  +  ri+l  +  ri+2  +  • ••  +  ri+n-l 
n 


Consequently,  the  next  rank  available  for  use  in  the  sequence  is  r,,_, 

i-rn 

In  either  case,  the  resultant  sequence  of  numbers  presents  the  numerical 
rank  of  each  value  in  the  original  time-series  relative  to  all  others  in 
that  series. 

Numerical  data  for  the  calculation  process  are  obtained  by  examining  the 

relationship  between  each  pair  of  ranks  in  the  modified  time-series  (r  , 

a 

r,  ,  r  ,  r.,  ...,  r  ),  beginning  with  the  element  r  .  The  u-1  pairs 

d  c  a  u  a 

containing  r  [(r  ,  r,  ),  (r  ,  r  ),  ....  (r  ,  r  )]  are  examined  first, 
a  a  b  a  c  a  u 

then  the  u-2  pairs  containing  r,  [eg,  (r,  ,  r.) . (r.  ,  r  )]  are 

b  Da  bu 

considered,  and  so  on  until  the  end  of  the  sequence  is  reached.  A  score 
of  +1  is  assigned  to  each  pair  having  the  ranks  in  ascending  order 
(eg,  the  pair  2,  6)  while  a  score  of  -1  is  assigned  to  each  pair  having 
the  ranks  in  inverse  order  (eg  the  pair  7,  5).  Scores  for  each 
individual  pair  are  summed  to  obtain  the  total  score  (S)  for  the 
modified  time-series.  If  there  are  no  tied  data  points,  the  Kendall 
rank  correlation  coefficient,  x,  is  defined  by  the  following: 


T  - 


_ total  score _ 

maximum  possible  score 


S 

^N(N-l)  ’ 
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where  N  -  the  number  of  measurements  in  the  time-series. 

When  there  are  tied  data  points,  the  Kendall  rank  correlation 
coefficient  is  defined  by: 


T 


_ s _ 

teNCM-l)]*5  teN(N-l)  -  T]*5 


where  7  •  k  £tj  (t^  -  1) , 


and 


the  number  of  elements  in  the  jth  group  of  tied  data  points. 


The  Kendall  rank  correlation  coefficient  can  have  any  value  from  +1  to 
-1.  If  the  series  of  ranks  is  a  monotonially  ascending  sequence  from  1 
to  N,  it  is  obvious  that  T  will  be  equal  to  1.  Similarly,  if  the  series 
of  ranks  is  a  monotonically  descending  sequence  from  N  to  1,  it  is  clear 
that  T  will  be  equal  to  -1.  The  value  of  T  observed  in  any  time-series 
may  be  used  directly  as  a  criterion  to  eliminate  highly  correlated  data 
sets.  For  example,  any  time-series  could  be  considered  unsuitable  for 
deconvolution  processing  if  T  >  T 


As  stated  at  the  outset,  the  Kendall  rank  correlation  coefficient  is 
used  to  detect  trends  in  time-series  data.  Consequently,  the  expected 
distribution  of  T  must  be  determined  for  samples  drawn  from  a  randomly- 
distributed  population.  Assuming  that  each  element  of  the  time-series 
has  an  equal  probability  of  occupying  any  position  in  the  temporal 
sequence,  then  the  expected  probability  distribution  of  the  parameter  T 
can  be  developed  for  every  sample  size  N.  Indeed,  such  tables  are 
available  -  even  for  moderately  large  sample  sizes  fie,  N  <  40).  For 
larger  samples  from  a  randomly-distributed  population,  the  parameter  T 
can  be  considered  as  a  normally  distributed  variate  with  a  mean  value  of 
zero  and  a  variance  of: 


ol  _  2(2N  +  5) 

9N(N-1)  ’ 


The  algorithm  employed  in  the  onboard  acoustic  data-processing  software 
calculates  the  quantity  ZTAU  (figure  5,  column  19)  which  expresses  the 
Kendall  rank  correlation  coefficient  actually  observed  (t)  in  terms  of 
the  parameters  for  the  expected  distribution: 


ZTAU  -  —  . 

T 

The  statistic  ZTAU  is  interpreted  and  used  in  much  the  same  manner  as 
the  terms  ZRUN  and  ZMSSD  described  previously.  Each  is  a  normally 
distributed  variate  with  a  mean  value  of  zero  and  a  unit  variance.  They 
are  used  in  conjunction  with  tabulated  values  for  the  normal 
distribution  to  confirm  or  reject  statistical  hypothesis  relating  the 
time-series  under  consideration  with  a  known  (normal)  probability 
distribution. 
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3.2.4  Power  difference 

One  of  the  simplest  and  most  powerful  tests  for  transients  in  the 
beam-noise  data  Is  the  difference  between  the  average  power  level  and 
the  geometric  mean  intensity  level  (PRDIF).  This  can  be  expressed  as 


PRDIF 


10  log 


10 


oyio) 


i«i 


i-l 


where  is  the  ith  beam  noise  level  in  this  beam-noise  time  series. 

Experience  indicates  that  under  normal  conditions  PRDIF  is  approximately 
2  to  2.5.  Greater  values  than  this  will  occur  when  high  level 
transients  are  present  in  the  data.  This  results  from  the 
characteristic  of  the  average  power  level  to  be  dominated  by  the  high 
level  components  while  the  geometric  mean  intensity  level  (dB  average) 
treats  the  data  as  an  arithmetic  sequence  in  decibels  and  is  relatively 
insensitive  to  a  few  high  level  data  points. 

3.2.5.  Spearman's  rank  correlation  matrix 

The  degree  of  association  between  two  series  of  measurements  can  be 

estimated  by  calculating  the  Spearman's  rank  correlation  coefficient, 

r  (ref. 7).  This  statistic  can  be  used  to  measure  the  correlation 
s 

between  noise-measurement  sequences  obtained  simultaneously  on  any 
two  beams.  Consequently,  if  any  two  series  of  beam-noise  measurements 
are  highly  correlated,  there  is  a  corresponding  probability  that  a 
common  noise-source  is  the  cause  of  both  sets  of  measurements.  An 
example  is  when  the  noise  from  the  towship  dominates  the  low-noise- level 
beams:  comparing  the  low-noise-level  beam  with  one  receiving  the  noise 
of  the  towship  gives  a  high  Spearman's  rank  correlation  coefficient 
which  is  evidence  that  the  two  are  correlated  and  a  potential  problem 
exists. 

Spearman's  method  of  analysis  is  similar  to  that  used  in  calculating 
Kendall's  rank  correlation  coefficient  because  it  deals  with  the  ranks 
of  the  observations  and  not  with  the  magnitudes  of  the  measurements 
themselves.  It  is  also  a  distribution-free  (non-parametric)  statistic 
and  does  not  depend  on  any  assumptions  concerning  the  probability 
distribution  of  the  values  in  either  time-series. 

To  calculate  Spearman's  rank  correlation  coefficient,  the  observations 
in  each  set  of  beam-noise  measurements  must  be  replaced  with  their 
ranking  relative  to  all  other  observations  in  that  set  (ie,  the 
measurements  from  beam  '  j  *  are  ranked  separately  from  those  of  beam 
'k').  The  procedure  used  to  rank  each  set  of  observations  is  identical 
to  that  described  previously  for  Kendall's  rank  correlation  coefficient 
and  will  not  be  repeated  here.  The  ranks  of  corresponding  positions  in 
the  two  temporal  sequences  are  then  compared  pair  by  pair  (ie,  the  rank 
of  the  first  measurement  in  the  time-series  for  beam  'j'  is  compared 
with  the  rank  of  the  first  measurement  in  the  time-series  for  beam  ' k ' ; 
this  process  is  repeated  for  each  of  the  N  positions  in  the  temporal 
sequences).  Spearman's  rank  correlation  coefficient  is  calculated  by 
the  following: 
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j-N 

1  '  N(N2-1)  \  dj 
j-1 


where  d 


j 


the  difference  between  the  ranks  of  corresponding  elements  In 
the  two  temporal  sequences. 


Spearman's  rank  correlation  coefficient  (rg)  can  have  any  value  from  +1 


to  -1.  If  the  two  series  are  identical,  it  is  obvious  that 


will 


always  be  zero  and  that  rg  will  equal  +1.  Similarly,  if  one  time-series 
is  a  mirror  image  of  the  other,  it  can  be  shown  that  rg  will  always  be 

equal  to  -1.  Spearman's  rank  correlation  coefficient  can  be  used 
directly  to  identify  data  sets  that  appear  to  be  strongly  correlated. 
For  example,  data  from  two  beams  that  are  widely  separated  in  space 

would  be  suspected  of  contamination  if  r  >  r  (determined  by  level 

s  s  max  J 


of  confidence  desired). 


The  coefficient  may  also  be  used  as  a  test  for  significance,  since  the 
expected  distribution  of  rg  is  known  to  be  symmetrical  around  the  value 

0,  to  approach  the  normal  curve  as  N  becomes  large,  and  to  be  truncated 
at  -1  and  +1.  Typically,  the  test  for  significance  is  applied  to  test 
the  null  hypothesis  that  the  rank  correlation  coefficient  in  the 
population  is  zero,  or  we  may  say  that  the  observations  in  the 
population  are  independent.  Tables  for  rg  are  available  for  small 

sample  sizes  (N  <  20) (ref .8,9) .  For  large  sample  sizes  the  sampling 
distribution  is  close  enough  to  normality  that  the  normal  area  table  may 
be  used  to  find  the  probabilities.  In  this  case  the  variance  of  r  is 
given  by! 


As  in  previous  tests,  the  observed  coefficient  (rg)  may  be  expressed  in 
terms  of  the  expected  distribution  (0,  a2)  by  the  following: 

ZSP  -  J-. 
r 


In  conjunction  with  tabulated  values  for  the  normal  distribution,  the 
teat  statistic  ZSP  may  be  used  to  confirm  or  reject  the  independence 
hypothesis  at  any  desired  level  of  significance. 

The  correlation  coefficients  and  their  corresponding  levels  of 
significance  (ZSP)  are  displayed  in  a  rectangular  matrix.  Figure  6  is 
an  example.  The  top  row  and  left-hand  column  are  beam  numbers  (the  last 
two  elements  in  each  are  for  hydrophones).  The  elements  above  the  main 
diagonal  are  the  correlation  coefficients  (xlOO,  for  display  purposes) 
and  those  on  and  below  are  the  confidence  levels  in  standard  deviations 
(xlO,  also  for  display  purposes).  To  facilitate  picking  out  possible 
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problem  areas,  all  significance  levels  less  than  three  standard 
deviations  (30  in  the  matrix)  have  been  replaced  by  zeroes.  The 
remaining  elements  indicate  that  the  hypothesis  of  independence  can  be 
rejected  with  a  confidence  level  of  greater  than  or  equal  to  99.73X. 

The  most  Important  part  of  the  matrix  is  what  is  printed  below  the  main 
diagonal.  A  non  zero  element  Indicates  correlation  in  a  pair  of  beams 
as  a  beam  and  hydrophone.  A  zero  below  the  main  diagonal  is  indicative 
of  non  correlation,  at  least  at  the  99.73J5  level.  Furthermore,  a  non 
zero  number  (ie  35)  below  the  main  diagonal  darkens  more  area  than  a 
single  zero.  Hence,  the  rank  correlation  matrix  can  be  used  as  grey 
scale  indicator  of  correlation.  The  darker  the  matrix  is  below  the  main 
diagonal  the  more  wide  spread  the  correlation  and  potential  array 
degradation.  This  test  used  in  conjunction  with  the  beam  noise  level 
plot  and  multiple  beam  levels  (figure  3)  is  an  excellent  indicator  of 
sidelobe  suppression  level  as  beam  noise  contamination  level  (see 
Section  3.2.6  on  virtual  beams). 

3.2.6  Virtual  beams 

At  frequencies  below  design  frequency,  virtual  beams  will  be  produced 
automatically  by  the  FFT  beamformer.  These  are  the  beams  that  are 
formed  by  time  delaying  (in  a  time-domain  beamformer)  or  phase  shifting 
(in  a  frequency-domain  beamformer  -  FFT  beamformer  for  example)  the 
signals  from  the  array  of  hydrophones  more  than  is  necessary  to  form  an 
endfire  beam.  These  are  the  beams  in  figure  3(b)  that  have  beam  numbers 
less  than  14  or  greater  than  50. 

Since  the  virtual  beams  do  not  have  main  lobes  in  real  acoustic  space, 
they  are  usually  Ignored.  However,  they  contain  energy  from  various 
incoherent  sources  of  self-noise  and  coherent  acoustic  energy  through 
sidelobes  in  the  acoustic  domain.  Simple  statistics  calculated  from  the 
outputs  of  these  virtual  beams  can  be  used  to  help  determine  the  levels 
and  sources  of  self  noise,  the  operating  condition  of  the  system,  and 
the  levels  of  sidelobe  suppression.  These  statistics  can  also  be  used 
as  a  signal  for  the  existence  of  problems  and  as  a  guide  for  locating 
the  problems. 

There  are  at  least  three  ways  in  which  the  virtual  beams  receive  energy. 

(1)  Virtual  beams,  like  real  beams,  have  sidelobes  that  extend  into 
real  space.  If  the  sidelobe  rejection  of  the  beamformer  is  poor, 
strong  sources  from  acoustic  space  can  "leak"  acoustic  energy  into 
the  virtual  beams. 

(2)  There  could  be  energy  propagation  in  the  array  at  a  speed  lower 
than  the  speed  of  sound  in  the  sea.  This  energy  would  appear  on  one 
virtual  beam. 

(3)  There  is  energy  on  the  hydrophones  that  is  not  coherent  from  ore 
hydrophone  to  another.  This  energy  may  be  of  non-acoustic  (flow 
noise),  electronic,  or  mechanical  (shocks,  vibrations)  origin.  These 
Incoherent  noises  are  spread  among  all  the  beams,  real  and  virtual. 
However,  they  are  most  easily  spotted  in  the  virtual  domain  because 
there  is  normally  less  energy  there  to  mask  them. 

The  energy  in  the  virtual  beams  can  therefore  be  invaluable  for  quality 
checking  and  "grooming"  the  towed-array  system. 

Additional  discussion  of  virtual  beams  is  available  in  reference  3. 
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4.  BEAM  NOISE  ANALYSIS  PRODUCTS 

Once  the  quality  of  the  measurements  has  been  evaluated  as  described  in 


Section 

analyse 

3,  the  appropriate  data-processing  products  of  Table  1  are  used  to 
the  ambient  noise.  This  is  done  under  six  headings: 

(1) 

Omnidirectional  noise  levels 

(2) 

Horizontal  directionality 

(3) 

High  resolution  spatial  spectrum 

(4) 

Temporal  statistics 

(5) 

Spatial  statistics 

(6) 

Array  heading  roses 

The  purposes  and  procedures  for  these  analyses 
Sections  4.1  through  4.5  respectively. 

are 

described  below  in 

4.1 

Omnidirectional  noise  levels 

Measuring  the  omnidirectional  level  of  underwater  ambient-noise  is  not  a 
simple,  straight-forward  task  when  the  measurement  tool  is  a  towed  array. 
This  is  because  the  towship  is  always  present  and  the  array  is  being  towed 
at  speeds  of  up  to  five  knots  when  measurements  are  taken.  The  proximity 
of  the  towship  contributes  significantly  to  the  measured  noise.  The  noise 
received  from  the  towship  can  be  as  much  as  69  dB  at  1.5  kHz  and  89  dB  at 
300  Hz  if  the  distance  between  the  towship  and  the  array  is  1000  m  (which 
is  typical).  The  ambient-noise,  on  the  other  hand,  can  be  more  than  10  dB 
below  these  levels. 

The  technique  used  to  estimate  the  omnidirectional  level  of  the 
ambient-noise  is  the  same  as  that  used  to  estimate  its  horizontal 
directionality,  which  is  described  in  Section  4.2  below.  In  fact  the 
omnidirectional  level  is  obtained  by  integrating  (or  summing)  the  levels  of 
noise  over  all  azimuths  in  the  directionality  pattern.  The  towship  noise 
has  very  little  Influence  on  the  directionality  pattern,  since  the  beams 
dominated  by  towship  noise  are  eliminated  in  the  noise-directionality 
calculations.  Hence  the  omnidirectional  level  is  an  average  value  for  the 
period  over  which  the  measurements  are  taken. 

The  output,  which  consists  of  a  single  value  for  each  analysis  frequency, 
is  printed  as  part  of  the  output  of  the  horizontal  directionality 
measurement,  as  shown  in  figure  8. 

4.2  Horizontal  directionality 

4.2.1  Purpose 

The  major  objective  of  the  measurements  made  during  the  directionality 
polygon  manoeuvres  is  to  provide  appropriate  data  to  resolve  ambiguities 
and  to  estimate  the  directionality  of  persistent  background  noise.  The 
concept  of  a  "persistent  background-noise  directionality"  is  applicable 
only  when  all  the  factors  that  Influence  the  ambient  noise  are 
reasonably  stable  over  long  periods  (6  hours  to  6  months)  and  when  the 
factors  are  generally  repeatable  so  as  to  yield  approximately  the  same 
directionality  pattern.  Seasonal  patterns  in  either 
acoustic-propagation  or  in  shipping  lanes  create  conditions  that  satisfy 
both  these  requirements:  their  contributions  to  the  ambient-noise  level 
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remain  reasonably  stable  during  a  given  season  and  repeat  themselves  for  < 

the  same  season  in  successive  years.  If  such  conditions  exist,  the  | 

final  criterion  is  the  relationship  of  these  factors  to  the  noise  level 

at  a  given  site.  For  example,  for  a  site  at  the  intersection  of  several  ■ 

shipping  lanes,  the  nearby  ships  may  appear  randomly  distributed  in 

range  and  azimuth  and  the  concept  of  a  persistent  background-noise 

directionality  may  not  be  applicable.  At  a  site  one  hundred  miles 

removed,  however,  the  shipping  patterns  may  appear  well  defined  and 

support  the  concept  of  persistent  background-noise  directionality. 

Because  the  same  general  directionality  pattern  can  be  expected  to  occur 
whenever  future  conditions  are  about  the  same,  the  concept  of  a 
persistent  background  noise  is  the  key  to  limiting  the  number  of 
measurements  required  to  characterize  the  ambient  noise  in  a  given  area. 

This  also  reduces  the  time  needed  to  acquire  a  data  base  for  the 

validation  of  ambient-noise  models. 

4.2.2  Method 

Horizontal  directionality  is  estimated  by  the  Wagstaff's  Iterative 
Technique  (WIT)  algor ithm(ref . 10) .  This  is  a  constrained  iterative 
restoration  algorithm  that  restores  some  of  the  high-order  spatial 

harmonics  filtered  out  by  the  low-pass  spatial  response  of  the  array.  | 

Its  present  version  contains  corrections  for  array  tilt,  three- 

dimensional  beam-response  patterns,  and  noise-field  vertical  arrival 

structure(ref . 11) .  The  technique  is  basically  as  follows.  WIT  uses  the 

median  beam-noise  values  from  all  legs  of  a  directionality  polygon  to 

"guess"  the  horizontal  directionality.  The  beam  patterns  of  the  array 

on  the  various  headings  of  the  polygon  are  then  convolved  with  this 

guess  to  calculate  median  beam-noise  levels.  These  calculated  levels 

are  compared  with  the  corresponding  measured  levels  and  the  "guess"  is 

appropriately  modified  along  azimuths  of  disagreement.  This  | 

convolution,  comparison,  and  modification  procedure  is  repeated  until 

the  agreement  is  within  predetermined  limits.  Such  an  iterative  I 

procedure  actually  deconvolves  the  array  beam  patterns  from  the 

data.  In  some  cases,  where  the  noise-field  is  reasonably  stationary,  I 

the  measured  and  calculated  beam  responses  are  so  close  that  it  could  be 

said  that  the  array  cannot  tell  the  difference  between  the  estimated  and 

the  actual  noise-fields.  When  data  from  only  one  heading  (leg)  are  used 

the  result  Is  a  deconvolution  of  the  data. 

The  noise  levels  of  the  beam  contaminated  by  the  towship  can  be  up  to 
20  dB  greater  than  beams  considered  free  of  contamination.  Even  so,  the 
ability  of  WIT  to  discriminate  against  or  eliminate  this  contamination 
is  sufficient  to  produce  good-quality  estimates  of  the  noise-field 
directionality. 

The  output  of  WIT  is  a  polar  plot  of  the  estimates  of  the  horizontal 
directionality  derived  from  the  50  percentile  (median)  noise  levels 
(figure  8).  This  is  done  for  each  analysis  frequency  of  the  data  set 
and  obtained  by  combining  the  median  beam-noise  levels  of  each  leg  of 
the  directionality  polygons.  The  result  is  usually  plotted  in  a  polar 
format  such  as  figure  8  and  is  often  referred  to  as  a  "noise  rose".  In 
addition,  as  explained  under  Section  3.1(h),  polar  plots  of  the 
50  percentile  (median)  noise  levels  on  all  beams  receiving  acoustic 
energy  (figure  7)  are  generated  for  the  time-series  of  each  leg  so  that 
anomalous  response  patterns  can  be  easily  spotted. 
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Figure  8.  A  data-processing  output  :  horizontal  directionality  plots  and 
omnidirectional  noise  levels 


4.3  High  resolution  spatial  spectra 

4.3.1  Purpose 

The  beam-noise  versus  beam  numbers  data  can  be  processed  by  the 
Wagstaf f-Berrou  Broadband  (WB2)  algorithm(ref . 12)  to  obtain  high 
resolution  ambiguous  spatial  spectra.  This  technique  estimates  levels 
and  relative  azimuths  much  better  than  can  be  done  from  the  original 
beam-noise  data.  This  is  illustrated  by  a  typical  example  in  figure  9. 

4.3.2  Method 

The  WB2  technique  is  a  high-resolution  algorithm  for  estimating  the 
noise-field  which  produced  the  beam-noise  data.  WB2  operates  on  the 
spectrum  analysed  beam-noise  data  for  one  time  period  rather  than  the 
hydrophone  data.  The  technique  starts  with  an  initial  guess  of  the 
spatial  spectrum  and  convolves  it  with  the  beam  patterns  in  cosine 
space.  The  differences  between  the  estimated  beam-noise  levels  and  the 
corresponding  measured  beam-noise  levels  are  used  to  appropriately 
modify  the  guessed  spectrum  within  the  coverage  areas  of  those  beams 
showing  the  disagreement.  The  resulting  new  spectral  estimate  is 
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convolved  with  the  bean  patterns  again  and  compared  with  the  measured 
beam-noise  levels.  This  process  of  modification  of  the  spectral 
estimate,  convolution,  and  comparison  with  the  measured  beam-noise  data 
is  repeated  until  the  agreement  is  within  predetermined  limits  or  a 
fixed  number  of  iterations  is  reached,  which  ever  comes  first. 

The  WB^  algorithm  can  be  considered  an  optimized  version  of  WIT  which 
takes  advantage  of  simplifications  in  testing  and  spectral  estimate 
modification.  These  simplif ications  are  made  possible  by  performing  the 
estimation  procedure  in  cosine  space  where  the  beams  are  all  the  same 
size  and  shape  and  by  estimating  the  ambiguous  noise-field  spatial 
spectrum  for  only  one  array  heading  and  one  time  period.  Like  WIT,  a 
key  to  the  success  of  WB2  is  that  the  estimation  process  is  performed  in 
the  log  transformed  space  (decibels)  rather  than  in  beam  power  space. 


Figure  9.  WB2  results  for  high  resolution  spatial  spectra 
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4.4  Temporal  statistics 
4.4.1  General 

The  minimum  number  of  samples  that  gives  statistical  significance  for 
one  time-series  is  considered  to  be  about  50.  During  the  directionality 
polygon  manoeuvres  data  are  acquired  during  15  min  time  periods  (about 
one  sample  every  15  s)  on  each  heading.  Thirty  minutes  are  usually 
required  for  turning  and  stabilizing  the  array  after  each  leg  of  the 
polygon.  Thus  it  is  just  possible  to  acquire  a  complete  set  of 
beam-noise  time-series  data  for  each  polygon  leg  in  about  45  min. 

However,  by  collecting  data  over  much  longer  periods  it  is  possible  to 
power-average  the  time-series  over  up  to  ten  consecutive  samples  to 
generate  new  time-series  of  down  to  one-tenth  of  the  original  size, 
thereby  giving  time-bandwidth  products  (TxBWs)  of  up  to  10.  This 
extends  the  applicability  of  the  results  to  those  sonar  systems 
operating  on  a  time-bandwidth  products  TxBW  of  greater  than  one. 

Thus,  although  the  data  collected  during  the  directionality  polygon 
manoeuvres  are  also  used  to  analyse  the  temporal  statistics,  that  is  not 
considered  to  be  their  primary  purpose.  Instead,  separate  "statistics 
tows"  are  conducted  especially  for  that  purpose.  These  tows  last  for  up 
to  12  h  on  a  constant  heading,  allowing  the  number  of  samples  in  their 
data  sets  to  reach  the  present  processing  limitation  of  250.  This 
allows  the  creation  of  analysis  products  with  TxBWs  of  up  to  10. 

The  analysis  products  generated  to  characterize  the  temporal  statistics 
of  the  ambient  noise  are: 

(a)  10,  25,  50  (median),  75  and  90  percentiles  of  noise  levels  at 
all  beams  and  two  hydrophones.  Single  listing  for  each  analysis 
frequency  in  each  time-period,  as  shown  in  the  10%  to  90%  columns  of 
figure  5. 

(b)  Geometric  mean  intensity  levels  (decibel  average)  on  all  beams 
and  two  hydrophones.  Single  listing  for  each  analysis  frequency  in 
each  time-period,  as  shown  in  the  AVG  column  in  figure  5. 

(c)  Average  noise-power  levels  on  all  beams  and  two  hydrophones. 
Single  listing  for  each  analysis  frequency  In  each  time-period,  as 
shown  in  the  AVGPR  column  in  figure  5. 

(d)  Standard  deviations  of  all  beams  and  two  hydrophones.  Single 
listing  for  each  analysis  frequency  in  each  time-period,  as  shown  in 
the  STDEV  column  in  figure  5. 

(e)  Skew  and  kurtosis  for  all  beams  and  two  hydrophones.  Single 
listing  for  each  analysis  frequency  in  each  time-period,  as  shown  in 
the  SKEW  and  KURT  columns  in  figure  5. 

(f)  Plots  of  the  cumulative  distribution  functions  of  noise  levels 
for  selected  beams  or  groups  of  beams.  Done  for  each  analysis 
frequency  for  each  time-period  and  for  the  combined  data  set  from 
each  directionality  polygon,  as  shown  in  figure  (10)  in 
Section  4.4.2. 
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(g)  Polar  plots  of  the  estimates  of  noise  horizontal  directionality 
derived  from  the  10,  25,  50  (median)  percentile  noise  levels  for 
beams.  Done  for  each  analysis  frequency  for  the  combined  data  set 
from  each  directionality  polygon,  as  shown  in  figure  12  in 
Section  4.4.3. 

The  uses  and  methods  of  calculating  products  (a)  to  (e) ,  which  are 
listings  of  values  for  each  analysis  frequency  in  each  time-series,  are 
self-explanatory.  Products  (f)  and  (g)  are  plotted  products  however, 
and  are  described  separately  in  Sections  4.4.2  and  4.4.3  below. 
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Figure  10.  A  data-processing  output  :  plots  of  the  cumulative  distribution 
functions  of  noise  levels  for  selected  beams  or  groups  of  beams 


-  25  - 


WSRL-TN-47/88 


4.4.2  Cumulative  distribution  functions 

Cumulative  distribution  functions  (CDF)  of  noise-levels,  of  the  type 
shown  in  figure  10,  are  generated  from  the  time-series  by  accumulating 
all  the  unaveraged  (TxBW-1)  noise-levels  for  a  single  selected  beam  or 
for  specified  groups  of  beams.  For  data  from  the  directionality 
polygons  this  is  done  for  each  analysis  frequency  for  each  time-period 
and  also  for  the  combined  data  set  obtained  by  combining  the 
time-periods  of  all  the  legs.  The  results  can  thereby  be  applied  to 
estimating  the  acoustic  performance  of  a  horizontal  array  on  an 
arbitrary  heading. 

The  time-series  of  the  long  statistical  tows  are  processed  in  the  same 
manner  as  the  time-series  of  the  directionality-polygon  legs  except 
that,  as  explained  in  Section  4.4.1  above,  samples  were  grouped  to  give 
T  x  BWs  of  up  to  10.  The  effect  of  changing  the  T  x  BW  must  be 
understood  in  order  to  interpret  the  results  correctly.  This  effect  can 
be  illustrated  by  producing  listed  statistical  products  (of  the  type 
shown  in  figure  5)  for  the  same  data  set  but  for  different  T  x  BWs.  As 
an  example,  figure  11  shows  the  results  of  averaging  the  power  of  199 
data  points  from  selected  aft  beams  at  500  Hz.  The  top  section  of  the 
figure  is  for  no  averaging  (T  x  BW«1);  the  second  is  for  averages  of  two 
consecutive  points  (T  x  BW-2) ;  the  next  is  for  averages  of  five 
(T  x  BW-5) ,  and  the  bottom  is  for  averages  of  ten  (T  x  BW-10). 

It  is  seen  from  figure  11  that  as  T  x  BW  Increases,  the  intensity  of  the 
noise  levels  Increase  by  about  7  dB  at  the  lowest  percentile,  increase 
by  about  1  dB  at  the  median,  and  slightly  decrease  at  the  highest 
percentile.  The  geometric  mean  intensity  levels  are  about  2.5  dB  higher 
for  a  T  x  BW  of  10  than  they  are  for  a  T  x  BW  of  one.  These  shifts  to 
higher  levels  with  increasing  T  x  BW  result  from  the  logarithmic 
relationship  between  power  and  level,  in  which  high  values  dominate  the 
level  of  the  average  power  more  than  they  do  the  average  of  the 
individual  levels. 

The  differences  in  average  noise-power  level  of  a  fraction  of  a  decibel 
between  lowest  and  highest  T  x  BWs  result  from  199  data  points  being 
used  for  the  series  In  the  upper  section,  198  for  the  second,  195  for 
the  third,  and  only  190  for  the  series  in  the  lowest  section.  The 
standard  deviations  decrease  by  about  3  dB  from  the  lowest  to  the 
highest  T  x  BW. 

The  remaining  parameters  also  show  some  effect.  The  skew  has  increased, 
the  kurtosis  has  decreased,  and  the  number-of-runs  statistic,  the  mean- 
square  successive  difference  statistic,  and  the  Kendall's  rank 
correlation  coefficient  have  all  decreased  with  increased  T  x  BW. 
little  confidence  can  be  placed  in  these  latter  results,  however, 
because  of  the  small  number  (19)  of  samples  with  T  x  BW»10. 

4.4.3  Horizontal  directionality  estimates  from  percentile  noise  levels 

Estimates  of  the  horizontal  directionality  of  the  noise  derived  from  the 
10,  25,  50  (median)  percentiles  for  all  azimuths  are  generated  as  polar 
plots,  as  shown  in  figure  12.  The  50  percentile  (median)  curve  is 
simply  the  standard  horizontal-directionality  curve  and  is  applicable  to 
any  array  whether  it  is  linear,  circular,  spherical,  random,  etc., 
provided  that  the  beams  receive  most  of  their  noise  power  within  about 
20°  of  the  horizontal.  This  pattern  could  simply  be  convolved  with  the 
beam  pattern  of  the  array  to  obtain  estimates  of  median  beam  levels. 
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Figure  11. 


The  effect  of  changing  the  time  x  bandwidth  product  in  analyring 
statistical  data 
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Figure  12.  A  data-processlng  output:  polar  plots  of  the  estimates  of 
horizontal  noise  directionality  (in  dB//uPa2/Hz. °)  derived 
from  the  10,  25,  50  (median)  percentile  noise  levels  for 
beams 


The  10  and  25  percentile  values,  on  the  other  hand,  cannot  be  used  in 
such  a  straighf orward  manner.  The  values  in  the  10  and  25  percentile 
curves  are  obtained  by  the  power  sum  of  approximately  equal  low-level 
noise  from  an  ambiguous  beam  pair  (the  conical  beam  of  the  line  array 
intersecting  two  different  azimuthal  sectors)  and,  hence,  are  valid  only 
for  beam  pairs,  A  single  unambiguous  beam  in  the  direction  of  the  low 
values  would  measure  an  entirely  different  cumulative  distribution 
function,  one  with  a  larger  spread.  This  is  not  the  case  in  the 
directions  of  the  high  levels,  however,  since  the  noise  is  generally 
received  by  a  single  beam  (one  of  the  two  azimuthal  sectors  intersected 
by  the  conical  beam  of  the  line  array).  Hence,  keeping  this  one 
limitation  in  mind  and  making  the  appropriate  modifications  in  the 
calculation,  the  10  and  25  percentile  values  can  be  estimated  for  any 
array  by  convolving  the  beam  pattern  of  the  array  with  the  corresponding 
curve. 
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4.5  Spatial  statistics 

The  azimuthal  anistropy,  or  spatial  granularity,  of  the  ambient  noise 
measured  by  an  array  can  be  expressed  quantitatively  by  plotting  its 
spatial  cumulative  distribution  function,  as  shown  in  figure  13.  These  are 
called  azimuthal  anistropy  cumulative  distribution  function  ( AACDFl  plots. 
AACDF  plots  are  produced  for  each  analysis  frequency  for  selected  time- 
periods  from  the  statistical  tows  and  for  the  data  set  obtained  by 
combining  the  time-periods  from  each  leg  of  a  directionality  polygon. 

The  AACDFs  preserve  the  spatial  variability  of  the  noise-field  that  the 
horizontal-directionality  estimates  average  out.  It  is  a  statistic  that 
characterizes  the  noise-field  by  emphasizing  both  the  spatial  and  the 
temporal  influence  of  the  noise-field  on  the  array.  The  AACDFs  are 
applicable  only  to  a  horizontal  line-array.  The  horizontal-directionality 
plots  (figure  8)  described  in  Section  3.2,  on  the  other  hand,  are 
relatively  system-independent  and  the  variability  in  time  is  averaged  out. 

The  AACDF  plots  convey  a  considerable  amount  of  information  about  the 
spatial  distributions  of  the  noise  sources.  This  information  is  valuable 
for  predicting  the  acoustic  performance  of  towed-array  systems  and  for 
assessing  the  correctness  of  the  spatial  distribution  and  received  noise 
levels  in  noise  models.  The  AACDF  plots  can  be  used  in  system-performance 
estimation  by  providing  spatial  distribution  functions  of  the  50  percentile 
(median)  noise  levels.  This  is  done  for  beamwidths  (or  sector  widths) 
ranging  from  0.5°  to  10°.  Each  point  on  an  AACDF  plot  corresponds  to  a 
beamwidth  B  (read  along  the  vertical  axis),  a  percentage  of  all  azimuth 
spaces  S  (read  along  the  horizontal  axis),  and  a  median  noise-level  on  the 
beam  L  (read  along  the  noise- level  curve) .  The  meaning  of  these  three 
numbers  is  that  a  line  array  with  a  beamwidth  of  B°  will  measure  median 
noise  levels  of  less  than  L  dB  on  S%  of  azimuth  space.  It  is  a  statistical 
measure  of  the  heights  and  frequencies  of  occurrence  of  the  high  levels  of 
the  noise-field,  of  the  depths  and  widths  of  the  low  regions,  and  of  the 
distribution  of  both  high  and  low  regions  in  space.  This  information  is 
presented  In  terms  of  a  parameter  -  the  beamwidth  -  that  makes  sense  for  a 
line  array. 

The  utility  of  the  AACDF  for  validating  noise-models  is,  perhaps,  not  as 
evident  as  it  is  for  estimating  the  performance  of  a  towed-array  system. 
One  of  the  first  things  usually  done  when  validating  a  model  Is  to  compare 
the  measured  omnidirectional  noise  level  with  the  one  calculated  by  the 
model  when  the  modelled  acoustic  and  noise-source  environment  is  as  close 
to  the  real  one  as  possible.  Next,  a  comparison  between  modelled  and 
measured  vertical-array  output  helps  to  determine  whether  the  physics  of 
acoustic-propagation  and  noise-source  radiation  are  correctly  modelled. 
Comparisons  of  horizontal  directionality  determine  the  correctness  of  the 
gross  spatial  distributions  of  the  noise-sources.  For  example,  if  a 
shipping  lane  is  on  one  side  of  a  site  but  modelled  on  the  opposite  side, 
the  modelled  omnidirectional  levels  and  vertical  directionality  can  be  in 
agreement  with  measured  values  but  the  horizontal  directionality  cannot. 

When  the  shipping  lane  contains  a  given  number  of  ships,  the  modelled 
horizontal  directionality  will  still  agree  with  the  measured  horizontal 
directionality  whether  the  model  has 

(a)  the  same  number  of  ships  with  the  correct  source  levels, 

(b)  ten  times  as  many  ships  in  the  same  general  direction  with  one- 

tenth  the  source  level,  or 


(c)  one-tenth  the  ships  with  ten  times  the  source  level. 
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Figure  13.  A  data-processing  output  :  plots  of  the  azimuthal  anistropy 
cumulative  distribution  function  (AACDF) 


The  AACDFs  for  modelled  and  measured  results,  however,  can  agree  only  if 
the  spatial  distribution  of  ships  is  correct.  This  is  because  the  AACDFs 
are  generated  from  the  outputs  of  beams  of  varying  widths,  some  that  are 
narrow  enough  to  fit  between  ships  and  others  that  are  not.  Only  when  the 
azimuthal  spacing  between  ships  in  the  model  is  statistically  equal  to  that 
which  existed  during  the  measurement  can  the  AACDF  plots  be  similar, 
ie,  have  the  same  spacing  and  slopes  of  curves;  and  only  when  the 
combination  of  the  average  source  level  and  the  propagation  loss  gives  the 
correct  received  level  can  the  modelled  AACDF  curves  have  the  same  levels 
as  those  obtained  from  measurements.  Hence  the  AACDF  is  a  powerful  tool 
for  validating  noise-models. 

4,6  Array  heading  rose 

Plots  that  give  the  S/N  gain  improvement  over  broadside  performance  can  be 
generated  from  noise  roses  such  as  the  ones  in  figure  8,  since  the  average 
or  median  array  beam-noise  for  any  array  heading  and  beam  steering  angle  in 
a  given  detection  sector  can  be  estimated  from  the  noise  roses.  Beam-noise 
estimates  are  obtained  by  convolving  the  beam  patterns  of  the  arrav  with 
the  noise  rose.  The  mean  beam-noise  levels  of  the  beams  with  steering 
angles  within  the  detection  sector  are  averaged,  and  the  resultant  is 
subtracted  from  the  average  noise  level  of  a  broadside  beam  in  that  sector. 
The  result  is  plotted  along  the  azimuth  of  the  array  heading  and  represents 
the  S/N  gain  improvement  over  broadside  performance  for  detection  of 
targets  in  the  detection  sector  for  the  array  heading  along  that  particular 
azimuth.  When  this  calculation  is  carried  out  for  all  possible  array 
headings  and  plotted  in  polar  form,  "array  heading”  plots  such  as  are 
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Figure  14.  A  date-processing  output  :  array  heading  roses  for  a  noise-field 
based  10  dB  up  in  the  northern  half  space 
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illustrated  in  figure  14  are  obtained.  Each  of  these  plots  is  for  a 
different  detection  quadrant.  The  specific  quadrant  for  each  is  denoted  by 
the  cross-hatched  area.  Sectors  with  other  orientations  and  widths  could 
have  been  chosen;  however,  these  sectors  (the  four  principal  quadrants)  are 
sufficient  to  illustrate  the  value  of  the  array  heading  roses. 

The  array  heading  roses  in  figure  14  correspond  to  a  noise-field  which  has 
10  dB  more  noise  in  the  northern  half-space  than  in  the  southern  half¬ 
space.  The  points  on  the  curves  outside  the  large  circles  (0  dB)  of  the 
array  heading  roses  correspond  to  array  headings  for  which  there  is  a  S/N 
gain  Improvement  over  broadside  performance  for  detection  of  targets  in  the 
detection  sector  (cross-hatched  sector) ,  while  the  associated  value 
indicates  how  many  decibels  of  improvement  to  expect.  The  points  on  the 
curves  inside  the  large  circles  (0  dB)  correspond  to  array  headings  for 
which  the  S/N  gain  performance  is  degraded  relative  to  broadside 
performance,  while  the  associated  value  gives  the  expected  degradation  in 
decibels.  By  this  means  these  plots  aid  in  the  selection  of  the  best  array 
headings  for  detection. 


5.  CONCLUSION 

This  report  has  described  the  details  and  usage  of  the  Towed  Array  Noise 
Analysis  Package  (TANAP).  Products  proven  useful  in  the  quality  assessment  of 
towed  array  data  have  been  shown  and  explained.  Once  the  quality  of  data  is 
assured,  the  same  package  can  be  used  to  analyse  the  combined  beam-noise  data 
to  obtain  estimates  of  the  horizontal  noise  directionality.  Finally  the  array 
heading  rose  is  shown  to  provide  assistance  in  using  a  towed  array  for 
acoustic  signal  detection  in  an  azlmuthally  dependent  noise  field. 

The  appendices  show  how  the  software  is  implemented  on  an  IBM  3033. 

The  package  presented  here  can  easily  be  enhanced  by  additional  products  and 
therefore  should  provide  a  solid  software  foundation  for  most  users  of  towed 
acoustic  array  systems. 
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APPENDIX  I 

IMPLEMENTATION  ON  AN  IBM  3033  COMPUTER 


1. 1  Overview 

The  purpose  of  these  Appendices  I  to  III  is  to  provide  basic  user 
documentation  for  the  beam-noise  data-processing  and  reduction  system 
implemented  on  the  IBM  3033  computer.  Included  herein  are  brief 
descriptions  of  each  program  and  subroutine. 

The  input  data  for  the  processing  program  are  written  on  disk  by  a 
beamformer.  The  disk  contains  the  beam  complex-coefficient  data  resulting 
from  the  double  FFT  beamforming.  The  software  for  this  task  is  called 
RAW. BEAMS. FORT  and  is  called  by  RAW. BEAMS . CNTL.  Data  are  calibrated  here. 

1.2  Program  and  subroutine  relationships 

Table  1.1  contains  an  overview  of  the  relationships  between  the  programs 
used.  The  ten  main  programs  call  subroutines  that  are  both  internal  and 
external.  The  table  attempts  to  Indicate  this  so  that  the  programmer 
seeking  a  particular  listing  will  know  where  to  look.  Note:  Some  internal 
subroutines  are  used  in  more  than  one  external  subroutine. 

1.3  Major  Noise  Program  Segments 

The  Towed-Array  Noise  Analysis  Packages  (TANAP)  consists  of  the  following 
major  segments: 

(a)  RDLMA 

RDLMA  selectively  reads  data  from  a  beamformer  disc  (BEAM)  and  creates 
data  output  for  4  frequencies.  The  data  are  formatted  in  a  data  file 
(NSTAT),  which  is  used  by  the  NSTAT  programs.  ICONT,  HYDAT,  PHASPL  are 
also  output. 

(b)  HYPLT 


Reads  HYDAT 

4  frequencies. 

and 

BEAM(headers) 

and 

plots 

hydrophone 

levels 

for 

(c)  PAPLT 

Reads  PHASPL 

4  frequencies. 

and 

BEAM(headers) 

and 

plots 

hydrophone 

phases 

for 

(d)  NSTAT 


NSTAT  processes  the  data  (NSTAT,  ICONT  BEAM(headers)  files  and  prints 
the  output  in  a  readable  format.  During  processing,  if  the  option  has 
been  selected,  the  program  prepares  a  data  file  (CDFST)  for  use  by  the 
Cumulative  Distributed  Function  program.  Output  to  files  PLDAT,  LONGS, 
F10,  F25,  NOISE  for  later  processing. 

(e)  BEAMPLT 

Reads  PLDAT  and  BEAM(headers)  and  plots  beam  levels  for  •'  frequencies. 

(f)  CDFPL 

CDFPL  plots  the  data  stored  In  the  data  file  (CDFST). 
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fg)  WB2 

Reads  PLDAT,  calculates  and  plots  high  resolution  spatial  spectra  using 
the  WB2  technique. 

(h)  NOISE 

NOISE  is  a  noise-ambiguity  resolution  program  that  determines  the 
ambient-noise  level  and  the  direction  from  which  the  noise  is  coming. 
The  input  data  for  this  program  are  prepared  by  the  NSTAT  program  and 
written  in  the  NOISE  file.  Additional  data  are  found  in  ICONT  and 
BFAM(headers) .  As  part  of  the  output  for  this  program  a  data  file 
(AMDAT)  is  prepared  for  use  by  the  NOIS7  program.  Files  L10  and  L25  are 
also  output. 

(i)  N0IS7 

N0IS7  is  a  noise-field  directionality  plot  program  that  plots  the 
results  of  the  NOISE  program.  It  uses  AMDAT  as  input. 

( j )  AHROZ 

AHROZ  uses  the  AMDAT  file  to  calculate  and  plot  array  heading  roses  for 
given  sectors  of  search. 

1.4  Files  Used 

The  following  files  are  used  in  TANAP. 


FILE 

UNIT 

FORMAT 

LONGS 

20 

BEAMS  FOR  LONG  STATS  RUN 

VB 

A4 

NSTAT 

21 

BEAM  OUTPUTS  X  FREO  X  ITERATIONS 

VB 

A4 

PLDAT 

22 

BEAM  DATA  FOR  PLOTTING  AND  WB2 

FB 

TEXT 

PHASPL 

23 

HYD  AVERAGE  PHASE 

VB 

A4 

ICONT 

24 

DATA  HEADINGS  -  USEFUL  CONSTANTS 

VB 

A4 

F10 

25 

10  ZILE  FOR  EACH  PERIOD/ALL-NOISE 

FB 

TEXT 

F25 

26 

25  ZILE  FOR  EACH  PERIOD/ALL-NOISE 

FB 

TEXT 

RIO 

27 

R25 

28 

L  10 

29 

10  PERCENTILES 

VP 

A4 

1.25 

30 

25  PERCENTILES 

VB 

A4 

AMDAT 

31 

PLOTTING  DATA  FOR  NOIS7 

VB 

A4 

NOISE 

32 

VARIABLES  AND  MEDIAN  BEAMS  FOR  WIT 

FB 

TEXT 

HYDAT 

33 

HYD  AVE.  MAG. 

VB 

A4 

CDFST 

34 

DATA  FOR  CDFPL 

VB 

A4 

CDFT1 

35 

DATA  FOR  CDFPL 

VB 

A4 

INPUT  BEAMS  10 
FROM  DJK 
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APPENDIX  II 

SIMPLIFIED  EXAMPLE  OF  A  RUN  PROCEDURE 
II. 1  Run  procedure 

The  beam-noise  data-processing  programs  can  be  run  interactively  on  the 
IBM  3033  computer.  The  following  example  illustrates  chis  and  can  be  used 
as  a  guide  if  the  operator  responses  are  modified  to  reflect  the  parameters 
that  characterize  the  type  of  processing  desired. 

In  the  following  run  example  the  instructions  and  answers  follow  each  other 
on  seperative  lines.  To  analyse  noise  data  the  following  interaction 
between  the  operator  and  computer  takes  place,  commencing  with  command 
TANAP: 
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SSI 

ENTER  S  DEPTHS  IN  METERS  (EX.  2tl,23S.2SI,2S5,2?l. >  *  - 

lB.ll.ll.ll.lt 

ENTER  LATITUDE  ( DECS, MINS, SEC )  i. 

-31 ,1,1 

ENTER  LONGITUDE  < DEGS, BINS, SEC )  i_ 

7 

151,1,1 


*«*  NADS  PARAMETERS  tn 

scu  -  I 

HEADINGS  •  t  1 

DEPTHS  -  ID  IS 

JULIAN  DAY  •  SSI 

LATITUDE  -  -31  DEG 

LONGITUDE  -  151  OCC 

DO  YOU  WISH  TO  USE  THESE  UALUES, 


•  AIN 
•  AIN 
V  OR  N? 


II 


II 


SCU  - 

HEADINGS  • 
DEPTHS  - 
JULIAN  DAY 
LATITUDE  • 
LONGITUDE  • 


1  1 

II  II 

SSI 

-31  DEG 
1S1  KG 


I  AIN 
I  AIN 


START  TINE 


SSI  13  8  SS 


NUNSER  Of  ITERATIONS  -  61 

ENTER  ARRAY  HEADING  (KG) 

7 

125 

ENTER  SIKLOK  LEUEL  IN  Dl,  EQ  31 

7 

SI 


ITA  • 

61 

IFLCMT 

1 

IFLCNT 

S 

IFLCNT 

3 

IFLCKT 

4 

IFICHT 

5 

IFLCNT 

SI 

IFLCNT 

SI 

IFLCNT 

SS 

IFLCNT 

S3 

IFLCNT 

54 

IFLCNT 

55 

IFLCNT 

56 

IFLCNT 

5? 

IFLCNT 

58 

IFLCNT 

59 

IFLCNT 

61 

IFLCNT 

61 

PROCESSING  SEQUENCE 

1 i  RDLAA  READ  RAW  BEANS  AND  SPECTRA.  OUTPUT  AUER AGE 
BEANS  AND  PHONES 

Si  MYPLT  PLOT  PHONE  LEVELS  AFTER  RDLAA 

3i  PAPLT  PLOT  PHONE  PHASES  AFTER  RDLAA 

4t  NSTAT  TENPORAL  STATISTICS.  WRITES  NOISE  FILE  AFTER  RDLAA 
Si  ICANPIT  PLOT  BEANS  AFTER  NSTAT 

6 i  CDFPL  PLOT  CUNULATIUC  DISTRIluTIOH  FUNCTIONS  AFTER  NSTAT 
7t  UIS  UBS  PROCESS  OF  BEANS  AFTER  NSTAT 
Nt  NO  PROCESSING  CONTINUE  ON 

NULL:  NEXT  STEP  <S> 

ENTER  PROCESSING  STEP 
3 

ENTER  PLOT  REQUEST  NUN BE*  (9  FOR  HELP) 


HELP  ON  THE  PLOT  REQUEST  NUMBER  i 

NN  •  NO  PLOTS 

1111  ■  ALL  PLOTS.  ALL  PHONES 

DIGITS  ENTERED  MAY  BE  EITHER  *  1 *  OR  ••• 

4  DIGITS  MUST  BE  ENTERED 

STARTING  FROM  THE  LEFT,  THE  DIGITS  REPRESENT i 
PRAUG 
MEDIAN 
A  VC 
STDEV 


EXAMPLE  I  THE  PLOT  REQUEST  NUMBER  1111 

URL  GENERATE  PLOTS  OF  T*  PHONE  DATA  FOR  AUGPR.AVG  AND  STKV 


ENTER  PLOT  REQUEST  NUMBER  (9  FOR  HELP) 
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START  Of  PLOT.  OPTIONS’ 


TOP  CURVES! 

AVOPR . 

MEDIAN - 

AVC . 

FOR  ABSOLUTE  LEVELS.  ADD  'X'  DB  TO  ABOVE 

BOTTOM  CUAVESi 

STDEV - 

SET!  2581 3855 


PHONE  NUMBER  PHONE  NUMBER 


FREQ  =  25.00 


FREQ  =  50.00 


X  =32.30  OB 


X  =35.20  DB 


PK0CESS1IW  SEQUENCE 


ti  ROUIA  KM  MU  KM)  AM  SPECTRA.  OUTPUT  AUERAK 
KM)  PM  PHONES 

21  HVPLT  PLOT  PHONE  LEVELS  PETE*  RDLNA 
31  PPPLT  PLOT  PHONE  PHPSES  PETER  RDLNA 

<1  MSTPT  TEMORAL  STATISTICS,  UNITES  NOISE  FILE  AFTER  POLNP 
SI  IEPNPLT  PLOT  KPNS  AFTER  NSTPT 

SI  COfPL  PLOT  CUNULPTIWE  SISTRISUTION  FUNCTIONS  Pf TIP  NSTPT 

ti  uta  usa  process  of  scans  after  nstat 

m  NO  PROCCSSINC  CONTINUE  ON 

HULL l  NEXT  STEP  (31 


ENTER  PROCESSINC  STEP 


5?8jP3MK?5?«8!:3 

n«ncv<r  i-oin 

«n  *«*«  «  —  **  •*  •« 

r»  ok  <ii  omv  ai  n  uma  m  •  «i  ♦  tf>r»  m  m 

SKWMamsssirjji 

r*  •«  m  •<  m  «  nm  *  ni  •  m  «•  u>  r- n*  «u 

mitUmiimii 

(v+u<suic*  men 

=2*it*?»Pi!x«sM55 

’*","’"'"'*"232522aS228SSIS58»!;SI8!<S»SSlli*P;R 
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APPENDIX  III 

PROGRAMS  AND  SUBROUTINES  USED  IN  THE  TOWED  ARRAY 
BEAM-NOISE  ANALYSIS  SOFTWARE 

111.1  Statistical  calculations,  programs  and  subroutines 

NSTAT  Program,  See  Appendix  1.3. 

R3LMA  Subroutine,  interactively  sets  up  parameters  for  NSTAT. 

Selectively  reads  data  in  NSTAT  file  for  NSTAT  program  to  print 
out  statistics  tables. 

RDLMA  Program,  See  Appendix  1.3. 

STAT2  Subroutine,  performs  statistical  calculations  on  all  beams 
including  the  virtual  beams. 

STAT5  Subroutine,  calculates  Spearman's  rank  correlation  coefficients. 

MEDAM  Subroutines,  calculates  median  value  of  input  data. 

MEDAN  Subroutine,  slightly  different  from  MEDAM. 

CLOZE  Subroutine  called  by  MEDAM  and  MEDAN,  finds  the  index  of  the 

element  closest  to  the  given  ranked  position. 

SZTAU  Subroutine,  calculates  Kendall's  Test  for  correlation. 

GROUP  Subroutine,  groups  data  and  computes  cumulative  distribution 

function. 

JRANK  Subroutine,  places  the  rank  at  an  integer  vector  in  another 

integer  vector.  The  original  vector  in  unchanged. 

STAT6  Subroutine,  computes  cumulative  distribution  function. 

PROB  Subroutine,  calculates  azimuthal  arisotropic  cumulative 

distribution  function. 

PROBL  Subroutine,  plots  azimuthal  anisotropy  cumulative  distribution 

function. 

SETUP  Subroutine  used  in  probability  plots. 

HYPLT  Program,  See  Appendix  1.3. 

BEAMPLT  Program,  plots  average  power,  dB  average,  median  standard 

deviation  and  power  differences  of  the  beam  noise. 

PAPLT  Program,  see  Appendix  1.3. 

CDFPL  Program,  see  Appendix  1.3. 

111. 2  WIT  related  programs  and  subroutines 

NOISE  Program,  driver  for  function  involving,  WIT.  "How  many 

frequencies"?  "How  many  legs"?  "Enter  number  of  bad  phones". 
Use  median  or  dB  average,  what  beams  to  skip,  AACDF's,  etc. 

N0IS1  Subroutine,  setup  data  and  run  parameters  for  the  WIT  algorithm. 
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N0IS2 

CONV 

MSOD 

NOIS5 

BWCAL 

SMART 

SMSEG 

NOIS6 

NOXS7 

FLDPL 

LEGPL 

III. 3 

WB2 

SMUZ 

CNVLV 

RRESP 

OREST 

NWDAT 

WITPL 

AHROZ 

AMRES 

BPLOP 

III. 4 

DDNAD 

BWTST 

BHDX 


Subroutine,  driver  for  WIT  algorithm. 

Subroutine,  performs  convolution  for  WIT. 

Subroutine,  performs  testing  and  field  modification  for  WIT 
algorithm. 

5“^0Utlnc>  calculates  bounds  (HIBAR  and  LOBAR)  used  in  subroutine 
MSQD. 

Subroutine,  modifies  beam  widths  and  beam  center  angles  for  use  by 
WIT.  Contains  correction  factor  for  less  than  64  hydrophones. 

Subroutine  to  calculate  SMART  results. 

Subroutine,  sets  up  data  and  calls  SMART  to  get  10  and 
25  percentile  directionalities. 

Subroutine,  writes  noise  field  from  WIT  and  beam  data  to  plot 
file.  v 

Program,  see  Appendix  1.3. 

Subroutine,  plots  the  noise  field  from  WIT  on  polar  axes. 

Subroutine,  plots  the  leg  beam  noise  data  on  polar  axes. 

WB2  related  programs  and  subroutines 

Program,  see  Appendix  1.3. 

Subroutine,  a  data  smoother. 

Subroutine,  does  convolution. 

Subroutine,  calculates  beam  responses. 

Subroutine,  calculates  first  estimate  of  noise  field. 

Subroutine,  reads  data  for  plotting  WB2  from  file  PLDAT. 

Subroutines,  plots  WB2  result. 

Program,  see  Appendix  1.3. 

Subroutine,  reads  data  from  AMDAT  for  AHROZ. 

Subroutine  plots  noise  and  array  heading  roses. 

Set  up  and  general  purpose  subroutines 

Interactively  requests  nonacoustic  data,  ie  SCU  gain,  array 
heading,  Julian  day,  array  depths,  latitude,  and  longitude  and 
prints  them  out. 

Requests  interactive  input  of  frequency  array  type,  array  heading, 
shading  and  sound  velocity. 

Calculates  the  beam  numbers  of  the  useable  beams  for  the  FFT 
beamformer8. 


-  43  - 


WSRL-TN-47/88 


R2LMA  Subroutine,  gets  hydrophones  data.  Unscrambles  the  beams. 

Calculates  median  hydrophone. 

WEIGT  Function,  computes  inverse  Hann  shading  coefficients  for 

40  hydrophones  and  a  64  point  FFT. 

RDBEM  Subroutine  to  read  beam  data  called  by  STAT5. 

ASINE  Function  calculates,  the  arc  sine  function  and  range  of  result 
from  90  to  +90  to  0  to  180° . 

BUST  Calculates  beam  widths  and  centres. 

FRAME  Subroutine  used  in  plotting. 

LINE  Subroutine,  draws  line,  used  in  all  plotting  subroutines. 

XYPLT  Subroutine,  draws  a  line. 
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